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Abstract 

<D We present a projective symmetry group (PSG) analysis of the spinless excitations of Z2 spin 

liquids on the kagome lattice. In the simplest case, vortices carrying Z2 magnetic flux ('visons') 
are shown to transform under the 48 element group GL(2, Z3). Alternative exchange couplings can 
also lead to a second case with visons transforming under 288 element group GL(2,Z3) x D3. We 
study the quantum phase transition in which visons condense into confining states with valence 
bond solid order. The critical field theories and confining states are classified using the vison PSGs. 
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I. INTRODUCTION 



Frustrated quantum magnets have a long history in condensed matter physics. Due 
to their relative simplicity in comparison to itinerant electron systems, they provide an 
ideal playground to study strongly correlated states of matter and the effects of competing 
ground states. One long standing goal is to identify realistic systems that realize spin 
liquid ground states, i.e. strongly correlated states of localized spins that don't break any 
symmetries. A promising candidate is the antiferromagnetic spin- 1/2 Heisenberg model on 
the kagome lattice and its realization in nature in the from of the mineral Herbertsmithite. 1 ' 2 
Until recently several theoretical works based on dimer model approaches, series expansions 
as well as numerical calculations suggested that the ground state of this model is not a 
spin liquid, but a valence bond solid 3-7 , i.e. a state that doesn't break the spin rotation 
symmetry, but instead breaks a lattice symmetry. However, a recent DMRG study of Yan et 
al. 8 has provided striking evidence for a spin liquid ground state for the S = 1/2 Heisenberg 
antiferromagnet on the kagome lattice. Yan et al. found a gap to all excitations, and it is 
plausible that their ground state realizes a Z 2 spin liquid. 9-16 

The results of Yan et al. also indicate the presence of proximate valence bond solid 
(VBS) states in which the space group symmetry of the kagome lattice is broken, and 
the fractionalized excitations of the spin liquid are confined into integer spin states. The 
confinement quantum phase transition should be accessible in extended models with further 
neighbor exchange interactions, and numerical studies of such transitions can serve as a 
valuable probe of characteristics of the spin liquid. 

This paper shall classify elementary vortex excitations of the Z 2 spin liquid, carrying Z 2 
magnetic flux 9,17,18 , often called 'visons', which are analogous to the Abrikosov vortices of 
BCS superconductors 19 (after electromagnetism is replaced by strong coupling to a compact 
U(l) gauge theory). We will compute their projective symmetry group 20 (PSG), and their 
spectrum using an effective frustrated Ising model. 21 Note that the Ising 'spin' has nothing 
to do with the S = 1/2 spin of the underlying antiferromagnet, and it is instead the creation 
or annihliation operator of the vortex excitation which is centered on sites of a lattice dual to 
that of the antiferromagnet. For the kagome antiferromagnet, the Ising model resides on the 
dice lattice, and the simplest effective Ising model has a degenerate momentum-independent 
spectrum. 4 ' 22 We shall show how the PSG constraints allow a systematic analysis of further 
neighbor interactions in the effective Ising model. Such extended interactions must gener- 
ically be present, 23 and they lead to well-defined vison states with a finite effective mass. 
Depending upon the values of these effective interactions, we find 2 possibilities for the vison 
states: the simplest case has them transforming under the 48 element group GL(2,Z 3 ), the 
group of 2 x 2 matrices with non-zero determinant whose matrix elements belong to the field 
Z 3 . The more complex case has visons states of the 288 element group GL(2,Z 3 ) x D 3 . 

Armed with this description of the vison states, will propose quantum field theories for 
the confinement transitions of Z 2 spin liquids on the kagome lattice. These transitions 
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FIG. 1: Visualization of two different valence bond solid states as dimer coverings of the kagome 
lattice. Each dimer represents a frustrated bond in the ordered phase of the corresponding Ising 
model on the dual dice lattice. The left pattern, representing the VBS If phase (see text), is a 
hardcore dimer covering with a 12-site unit cell that maximizes the number of perfectly flippable 
diamonds (highlighted in gray). Note however, that the choice of diamonds is not unique because 
the dimer pattern is symmetric with respect to 7r/3-rotations around the hexagons marked by a 
black dot. The pattern on the right represents the VBS 1a phase (see text), has a 36-site unit cell 
and maximizes the number of perfectly flippable hexagons (highlighted in gray). 

are associated with condensation of visons, 18,21 ' 23 ~ 27 and are expressed in terms of a multi- 
component 'relativistic' scalar field; the field theory with GL(2,Za) symmetry appears in 
Eq. (3.8). These field theories also place constraints on the specific patterns of spatial broken 
symmetry in the confining VBS states found next to the quantum critical point and we will 
present phase diagrams illustrating these states. Visualizations of two possible VBS states 
are shown in Fig. 1. The left VBS pattern in Fig. 1 appears for the simplest case of the 
GL(2,Z3) visons, and it is interesting that it is closely related to the "diamond pattern" 
which is enhanced in the numerical studies of Yan et al. 8 . The right pattern in Fig. 1 is one 
of many possible VBS states for the GL(2, Z 3 ) x D 3 visons, and maps to the "honeycomb" 
VBS states found in earlier studies 3-7 . 

We will begin in Section II with a review of the basic characteristics of Z2 spin liquids and 
of their vison excitations. A key property of a vison is that it picks up a Aharanov-Bohm 
phase of ir upon encircling every 5 = 1/2 spin on the sites of the antiferromagnet; 18 ' 21 ' 24 ' 25 , 
as will be described in Section II. 

Section III contains our main new results. We begin with general effective theories of 
vison motion which incorporate the Aharanov-Bohm phase 23,26 of tt: for the kagome anti- 
ferromagnet, these are conveniently expressed in terms of an effective, fully-frustrated Ising 
model on the dice lattice. A key parameter in this Ising model is the sign of a particular 
next-nearest-neighbor interaction. A 'ferromagnetic' sign leads to the simpler vison PSG, 
and is discussed in Section III A; the more complicated 'antiferromagnetic' case is discussed 
in Section III B. In both cases, we use the vison PSG to present quantum field theories for 
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the confinement transitions, and discuss renormalization group analyses of these quantum 
critical points. The field theories are also used to classify the patterns of lattice symmetry 
breaking in the confining valence bond solid states. 

In Appendix A the explicit calculation of a visons Berry phase can be found. The results 
of a PSG analysis of visons for different lattice geometries are summarized in Appendix B. 

II. Z 2 SPIN LIQUIDS AND VISONS 

We begin by a review of the basic properties of Z 2 spin liquids, following the description 
in Ref. 11. 

It is convenient to describe the S — 1/2 spins Si using Schwinger bosons 28 bi a (a =t,i) 

1 . 

S i = ifiJa&Pi (2.1) 

where a are the Pauli matrices, and the bosons obey the local constraint 

S>U« = 1 (2-2) 

a 

on every site i. Our analysis below can be easily extended to gapped Z 2 spin liquids obtained 
from the Schwinger fermion formulation, 10 > 14 ~ 16 but we will only consider the Schwinger boson 
case for brevity. 

The Z 2 spin liquid is described by an effective boson Hamiltonian 

Hb = ~ E Qa £ «$A + H.c. + A J2 b lbia, (2.3) 

i<j i 

where e is the antisymmetric unit tensor, A is chosen to satisfy the constraint in Eq. (2.2) 
on average, and the Qij = —Qji are a set of variational parameters chosen to optimize the 
energy of the spin liquid state. Generally, the Qij are chosen to be non-zero only on near 
neighbor links. The £ Z 2 ' character of the spin liquid requires that the links with non-zero Qij 
can form closed loops with an odd number of links, 9-11 . If the state preserves time-reversal 
symmetry, all the Qij can be chosen real. 

This effective Hamiltonian also motivates a wavefunction of the spin liquid 13,29 

|SL> = V exp fa ^b\ a b))j |0>, (2.4) 

where |0) is the boson vaccum, V is a projection operator which selects only states which 
obey Eq. (2.2), and the boson pair wavefunction = —fji is determined by diagonalizing 
Eq. (2.3) by a Bogoliubov transformation (see Eq. (A5) in Appendix A). 
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FIG. 2: A vison on the kagome lattice. The center of the vison is marked by the X. We have 
sgn(Q^) = — sgn(Qij) only on the links marked by the wavy lines. 

The Schwinger boson approach also allows a description of the vison excited states. 9,19 
We choose the vison state as the ground state of a Hamiltonian, obtained from % by 
mapping — > (see Eq. (Al)); then the vison state \^f v ) has a wavefunction as in 
Eq. (2.4), but with — > (see Eq. (A4)). Far from the center of the vison, we have 
\Qij\ = \Qij\i while closer to the center there are differences in the magnitudes. However, 
the key difference is in the signs of the link variables, as illustrated in Fig. 2: there is a 
'branch-cut' emerging from the vison core along which sgn(Q^) = — sgn(Qjj). This branch- 
cut ensures that the Z 2 flux equals -1 on all loops which encircle the vison core, while other 
loops do not have non-trivial Z 2 flux. Such a configuration of the is expected to be 
a metastable solution of the Schwinger boson mean-field equations, representing a vison 
excitation. 

The previous solution for the vison state 19 was obtained using a continuum field theoretic 
representation of the Schwinger boson theory, valid in the limit of a small energy gap towards 
spinful excitations. In principle, such a solution can also be obtained by a complete solution 
of the Schwinger boson equations on the lattice, but this requires considerable numerical 
effort. Here, we illustrate the solution obtained in the of a large spin gap. In the large 
spin gap limit, 30 we can integrate out the Schwinger bosons, and write the energy as a 
local functional of the Qij. This functional is strongly constrained by gauge-invariance: for 
time- independent Qij, this functional takes the form 

E[{Q l3 }] = -Y J Lm 2 + ^\Q l A+K £ QaQ%--Qu (2-5) 

i<j ' even loops 

Here a, /3, and K are coupling constants determined by the parameters in the Hamiltonian 
of the antiferromagnet. We have shown them to be site-independent, because we have only 
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FIG. 3: A vison on the triangular lattice (a similar solution is expected on the kagome lattice). 
The center of the vison is marked by the X. The wavy line is the 'branch-cut' where we have 
sgn(Q^) = — sgn(Qjj) only on the links crossed by the line. Plotted is the minimization result of 
with a = = —2,K = —0.5. Minimization is done with the cluster embedded in a 
vison- free lattice with all nearest neighbor links Qij. The numbers are (Qij — Qij) arid the thickness 
of the links are proportional to (QV. - Qij) 1 / 2 . K < gives rise to the zero flux state while K > 
favors the tt flux state where the unit cell is doubled. The zero and it flux states without the vison 
have been studied previously. 11,32 

displayed terms in which all links/loops are equivalent; they can depend upon links/loops 
for longer range couplings provided the full lattice symmetry is preserved. We describe the 
results of a minimization of on the simpler case of the triangular lattice in Fig. 3. 

The magnitudes of are suppressed close to the vison, and converge to Qij as we move 
away from the vison (modulo the sign change associated with the branch cut), analogous to 
the Abrikisov vortices. Despite the branchcut breaking the 3-fold rotation symmetry, the 
gauge-invariant fluxes of Q\j preserve the rotation symmetry. 

Let us now consider the motion of a single vison. A key ingredient is the Berry phase a 
vison accumulates while moving through the background spin liquid. The gauge-invariant 
Berry phases are those associated with a periodic motion, and so let us consider the motion 
of a vison along a general closed loop C. We illustrate the simple case where C encloses 
a single site of the triangular lattice antiferromagnet in Fig. 4. The Berry phase for this 
periodic motion can be computed in a manner analogous to that presented in Section III. A 
of Ref. 31 for a monopole in a U(l) spin liquid; details appear in Appendix A. The final 
gauge-invariant Berry phase turns out to be given by the gauge-transformation required to 
map the final state to the initial state. The analysis in Fig. 4 shows that the required gauge 
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FIG. 4: Periodic motion of a vison around a closed loop C on the triangular lattice. Here C encloses 
the single site marked by the filled circle. The wavy lines indicate sgn(Qy) = — sgn(Qy), as in 
Fig. 2. The bottom state is gauge-equivalent to the top state, after the gauge transformation 
bia —bia only for the site i marked by the filled circle. 

transformation is 

bi a — > — h a , for i inside C 
bia b ia , for i outside C. (2.6) 

By Eq. (2.2), each site has one boson, and so the total Berry phase accumulated by \^ v ) is 

7r x (number of sites enclosed by C) , (2.7) 

as recognized in earlier works 18 ' 21 ' 23 ' 24 ' 26 . It is also clear that for a spin S antiferromagnet, 
the Berry phase would be multiplied by a factor of 25. 

Finally, let us also mention the spinon states, although these will not play a role in the 
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subsequent analysis of the present paper. These are created by applying the Bogoliubov 
quasiparticle operator 7* (Appendix A) on the spin-liquid ground state; in this manner we 
obtain the spinon state 



where U ifJj is a Bogoliubov rotation matrix defined in Appendix A. We can now also consider a 
spinon well-separated from a vison, and describe the motion of a vison along a large contour 
C which encircles the spinon. First, consider the motion of a vison around the spinon 
state \£a) localized on the site t. Proceeding with the argument as above, the projection 
onto states which obey Eq. (2.2) now implies that the Berry phase for such a process is 
71 x ((number of sites enclosed by C) — 1). The transformation to the spinon state \[ia) 
will not change the result, provided the 'wavefunction' {U^ 1 *)^ is localized well within the 
contour C. Thus relative to the Berry phase in the case without a spinon, the vison acquires 
an additional phase of n upon encircling a spinon i.e. the spinons and visons are relative 
semions. 17 

III. PSG ANALYSIS OF THE FULLY FRUSTRATED ISING MODEL ON THE 
DICE LATTICE 

The analysis of Section II suggests a simple effective model for the vison fluctuations 
about the Z 2 spin liquid ground state. In the framework of the path-integral formations of 
Ref. 11, the visons are saddle-points of the Qij with 7r-flux, as shown in Fig. 3. Each saddle- 
point traces a world-line in spacetime, representing the time evolution of the vison. The 
Berry phase computation in Section II shows that this world-line picks up 7r-flux each time it 
encircles a kagome lattice site. We know further that two such worldlines can annihilate each 
other, because 2ir magnetic flux is equivalent to zero flux. So if we interpret the wordlines 
as the trajectory of a particle, that particle must be its own anti-particle, and has a real 
field operator. In this manner, we see that the fluctuations represented as the sum over all 
vison worldlines is precisely that of a frustrated Ising model in a transverse field 18,21,24 : 
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(2.8) 




(3.1) 



i<j 



where the product of bonds around each elementary plaquette is negative 



(3.2) 



plaq. 
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FIG. 5: Left: unit cell of the dice lattice consisting of the 3 independent sites labeled A, B and 
C. The A sites are 6-coordinated, whereas the B and C sites are 3-coordinated. The basis vectors 
u = (3/2, \/3/2) and v = (3/2, — \/3/2) are indicated by dashed arrows. Right: gauge choice for 
the fully frustrated Ising model on the dice lattice with a 12-site unit cell. Red thick bonds are 
frustrated, i.e. sgn(Jjj) = — 1. 

Also, we have not displayed the transverse-field term, because most of our symmetry con- 
siderations are restricted to time-independent, static configurations. 

In the following we will use a soft-spin formulation where the <f>j's take real values. This 
model is invariant under Z 2 gauge transformations <\>j — > J%j — > a^ajJij with a,j = ±1. 
For the case of the Z 2 spin liquid on a kagome lattice the dual Ising model lives on the dice 
lattice, shown in Fig. 5. The dice lattice has three independent sites per hexagonal unit 
cell, two of them are three-coordinated and one is six-coordinated. Flipping a dual Ising 
spin changes the flux through a plaquette on the kagome lattice by it, thereby creating or 
annihilating a vison. The magnetically disordered phase of the Ising model corresponds to 
a Z 2 spin liquid with deconfined spinon excitations, whereas the ordered phases describe 
different valence bond solids where the visons are condensed and fractional excitations are 
confined. 

In the following we study the confinement transitions of the spin liquid by construct- 
ing a Ginzburg-Landau functional that is consistent with the projective symmetry group 
(PSG), i.e. the combination of lattice-symmetry- and Z 2 gauge-transformations that leave 
the Hamiltonian (3.1) together with (3.2) invariant. In order to determine the PSG trans- 
formations we fix the gauge of nearest neighbor interactions as shown in Fig. 5, thereby 
obtaining a unit cell with twelve sites. The generators of the dice lattice symmetry group 
are translations by one of the two basis vectors, e.g. u, reflections about one axis, e.g. 
the x-axis, and tt/3 rotations about the central 6-coordinated site. Since our gauge choice 
is already invariant under rotations, we only need to determine the gauge-transformations 
corresponding to translations and reflections to specify the PSG. These are shown in Fig. 6. 

The dispersion relations of the soft-spin modes can be obtained directly from the Hamil- 
tonian (3.1). The corresponding action takes the form 




(3.3) 
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FIG. 6: Gauge transformations associated with translations by u and reflections about the x-axis. 
Shown is the unit cell of the fully frustrated dice model after a lattice translation by u (left) and 
after a reflection about the x-axis (right). The corresponding gauge-transformations consist of spin 
flips on the sites marked by blue points, which restore the original gauge pattern as in Fig. 5. 



where the summation over the sublattice indices i,j = 1...12 is implicit, J q denotes the 
Fourier transform of the interaction matrix J^-, and we have included a kinetic energy with 
frequency Q (which descends from the transverse field), and a mass term, m. As noted 
earlier 4 ' 22 , the frustrated Ising model on the dice lattice with nearest neighbor interactions 
gives rise to three flat bands (each being four-fold degenerate in our gauge choice with a 
12-site unit cell), which would result in infinitely many critical modes. In order to lift this 
degeneracy we include further interactions beyond nearest neighbors 23 that are consistent 
with the PSG. Physically, these interactions correspond to the hopping of visons beyond 
nearest neighbors. Different masses on the 3- and 6-coordinated sites would be allowed by 
the PSG, but they don't give rise to a momentum dependence of the modes and thus don't 
change the picture qualitatively. Out of the five possible additional interactions up to a 
distance of two times the nearest neighbor bond length only one is consistent with the PSG. 
This interaction, shown in Fig. 7, connects different 3-coordinated lattice sites and gives rise 
to a non-flat dispersion. The Fourier transform of the interaction matrix Jij takes the form 
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(3.4) 



Here we have abbreviated k x = t(l + e i2qu + e i2q ( u ~ v )), k 2 = t(l + e i2qu + e i2qv ) and k 3 = 

£( 1 + e i2qv + e i2q(v-u)) ) wh j ch 

are the additional terms coming from the next nearest neighbor 
interaction that is allowed by the PSG and t denotes the relative strength of the interaction 
with respect the the nearest neighbor coupling J = 1. Diagonalizing results in three 
dispersing bands, each being four-fold degenerate. Depending on the sign of the next-nearest 
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FIG. 7: Additional next-nearest neighbor interactions between 3-coordinated sites on the dice 
lattice (shown as dashed blue lines) that are allowed by the PSG. For illustrative clarity not all 
additional bonds in the unit cell are shown. All other bonds can be obtained from the ones that 
are shown via translations by the basis vectors u and/or v. 



neighbor interaction t, we get different transitions to valence bond solids that we are going 
to discuss in the following. 

A. Ferromagnetic n.n.n. interactions 

For ferromagnetic next-nearest neighbor interactions (t > 0) the dispersion minimum 
of the lowest band is at zero momentum q = 0. The corresponding eigenvalue of the 
interaction matrix is given by A + = (3t + V24 + 9t 2 )/2 and one choice for the four 
degenerate eigenvectors is 



,( 2 ) 



,(3) 



1 A + -l 



XT 1 A 6 A_ 







-1 A+ 

A_ 6 



o a^ o o £ =f 
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A+ A + A + 
6 6 6 



6/(A^ + 6) 
6/M + 6) 



with A_ = (3t — y/24: + 9t 2 )/2. This set of eigenvectors forms an orthonormal basis for the 
four critical modes at the transition to the confined phase. In the magnetically ordered 
state the magnetization <f)j(R) at lattice site R = 2mi + 2mv (n,m G N) and sublattice site 
j E {1, 12} is given by 

0,(11) = ]T ij n vf (3.5) 

n=1...4 

and is independent of the lattice site R for ferromagnetic n.n.n. interactions, i.e. the ordered 
VBS phases have a 12-site unit cell. The values of the four mode amplitudes ip n are obtained 
by minimizing the Ginzburg-Landau functional, which in turn is given by all homogeneous 
polynomials in the mode amplitudes ip n that are invariant under the PSG transformations. 
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In order construct these polynomials we have to determine how the mode amplitudes trans- 
form under the PSG. This can be done by looking at the transformation properties of the 
eigenvectors v^ n \ For example, the transformation properties of the mode amplitudes with 
respect to translations T u are determined via 

T u 0, = ^ T u vf = ^(T u ) mn ^ m) = £(T u ^) m i;< m) . (3.6) 

n n,m m 

The resulting PSG transformation matrices for the amplitudes of the four critical modes 
with respect to translations (T u ), reflections (I x ) and rotations (Rq) are given by 
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(3.7) 



We used the GAP program 33 to show that these three matrices generate a finite, 48 element 
subgroup of 0(4) which is isomorphic to GL(2,Z 3 ). We also determined that this group is 
isomorphic to the group x Cq\ in the classification of Conway and Smith. 34 

Next we determined the most general Ginzburg-Landau (GL) functional of the ip n . It 
turns out that there are three fourth order polynomials that are invariant under this group, 
thus our functional for the four mode amplitudes depends on the coupling constants, r, u, 
a and b, and takes the form 

C = (Wn) 2 + (dM 2 + rl£ + u^ n ) + a £ ^ 

n=l..A n<m 

+b [^1(^2^3 - ^2^4 + ^3^4) + ^2(^1^3 + ^1^4 - ^3^4) + ^3(^1^2 - ^1^4 + ^2^4) 
-fc + Wa + M)]. (3.8) 

At b = and a = 2u, this is just the well-known 4 field theory with 0(4) symmetry, whose 
critical properties are described by the extensively studied Wilson-Fisher fixed point. The 
b coupling breaks the 0(4) symmetry down to GL(2,Z 3 ). Remarkably, the renormalization 
group (RG) properties of just such a quartic coupling have been studied earlier by Toledano 
et al.; 35 they denoted this symmetry class as [D 3 /C 2 ; OfD 2 ], following the analysis of Du 
Val. 36 Toledano et al. found that the 0(4) fixed point was unstable, but were unable to find a 
stable critical fixed point at two-loop order. It would be interesting to extend the analysis of 
Eq. (3.8) to higher loop order, and search for a suitable fixed point by the methods reviewed 
in Ref. 37. 

Despite the difficulty in finding a suitable critical fixed point, we can assume the existence 
of a second-order quantum phase transition, and use Eq. (3.8) to determine the structures 
of possible confining phases. Minimizing this functional for the magnetically ordered phase 
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FIG. 8: Phase diagram of Eqn. (3.8) as a function of the two couplings a and b (here we have set 
u = 1). The VBS If phase is not reflection symmetric, whereas the VBS 2p phase is reflection 
symmetric. In the VBS 2p phase there is a crossover (indicated by the dashed line) from a phase 
where one of the ipnS is zero (left of the dashed line) to a phase where three of the tpnS are zero 
(right of the dashed line). 

r < gives rise to two possible phases, depending on the values of the two parameters a 
and b in the GL-functional (3.8). The phase diagram is shown in Fig. 8. 

In the VBS If phase the GL-functional (3.8) has 16 degenerate minima corresponding to 
different magnetization patterns. Note, however, that the magnetization itself is not a gauge 
invariant quantity. In fact, all degenerate minima give rise to the same bond patterns, which 
are related by simple lattice symmetry transformations. The bond pattern in the VBS If 
phase is symmetric under rotations by tt/3, translations by 2u and 2v but it is not reflection 
symmetric. The 16 minima thus correspond to 8 different bond patterns that are related 
by translations by u and v as well as by a reflection about the x-axis (the factor of two 
arises from a global spin flip symmetry). In Fig. 9 we show the bond pattern in the VBS If 
phase, i.e. we plot the gauge-invariant value of the bond-strength Jij&ifrj for every nearest 
neighbor bond. In this figure we mapped the bond pattern from the dice lattice back to the 
kagome lattice by assigning the value of the bond-strength on the dice lattice to the bond 
on the kagome lattice that intersects the dice lattice bond. Deep in the confined phase the 
bonds plotted in Fig. 9 thus represent the expectation value (cr^) of the Z 2 gauge field on 
the kagome bonds in the transverse field direction. 

One way to relate the symmetries of valence bond solid phases on the kagome lattice to the 
ordered phases of the frustrated Ising model is to make use of the mapping between hardcore 
dimer models and frustrated Ising models 38 , i.e. by putting a dimer on every bond of the 
kagome lattice that intersects a frustrated bond on the dual dice lattice. The dimer covering 
that is obtained in this way for the VBS If phase is shown in Fig. 1. In comparison to 
previously obtained dimer coverings of the kagome lattice 4,27 which maximize the number of 
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FIG. 9: Bond pattern in the VBS If phase. Plotted is the gauge invariant bond-strength Jijfiitfij 
for nearest neighbor bonds on the dice- lattice at a = b = 1, which has been assigned to each 
respective kagome bond. Black dashed lines indicate satisfied bonds (— Jy0j^>j < 0), red solid 
lines are frustrated bonds (—Jijcftifij > 0). The thickness of the bonds is proportional to the 
bond-strength. 

perfectly flippable hexagons on a 36-site unit cell, this dimer covering maximizes the number 
of perfect flippable diamonds on a 12-site unit cell. Note that this bond pattern is similar 
to the diamond pattern observed by Yan et al. 8 . It is important to note, that this mapping 
between the frustrated Ising model (3.1) and the corresponding hardcore dimer model is 
strictly valid only in the limit where the transverse field vanishes, i. e. deep in the ordered 
phase. In this regime our GL-approach is quantitatively not reliable, however. It is thus no 
surprise, that not all of our dimer coverings are hardcore coverings (see e.g. the right dimer 
covering shown in Fig. 1). The GL-approach works well to determine broken symmetries 
of VBS states as well as the critical properties close to the confinement transition, where 
it is sufficient to keep only the relevant terms up to fourth order in the functional. On 
general grounds it would be necessary to include all higher order terms in the GL-functional 
in order to obtain reliable hardcore dimer coverings in the limit of small transverse fields. 
Apparently, the VBS If phase is an exception to this rule. 

We note that the VBS lp state has recently been identified as the ground state of the 
deformed kagome lattice spin-1/2 antiferromagnet Rb2Cu3SnF 12 . 39 It is possible that the 
spin physics being discussed here played a role in the lattice distortion observed in this 
experiment. 40 Also in Zn-paratacamite, Zn^Cu^^OH^C^, there is a transition 41 between 
a spin-liquid phase near x — 1 to a distorted kagome lattice near x = 0, and it has been 
argued 42 that a 'pinwheel' VBS state, which is identical in symmetry to our VBS lp state, 
plays a role in the latter case. 
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FIG. 10: Bond patterns in the VBS 2p phase. Plotted is the gauge invariant bond-strength J^^j 
for nearest neighbor bonds on the dice lattice, shown on the corresponding kagome bonds. Black 
dashed lines represent satisfied bonds {—Jij4>i4>j < 0). In this phase there are no frustrated bonds 
at all. The thickness of the bonds is proportional to the bond-strength. 



In the VBS 2 F phase the GL-functional (3.8) has 8 degenerate minima and the corre- 
sponding bond-patterns have an additional reflection symmetry as compared to the VBS If 
phase. Moreover, there is a crossover from a phase where one of the ^ n 's is zero to a phase 
where three of the ifj^s are zero, indicated by the dashed line in Fig. 8. The corresponding 
bond patterns are shown in Fig. 10. 



B. Antiferromagnetic n.n.n. interactions 



If the additional next-nearest neighbor interactions are antiferromagnetic (t < 0) and 
smaller than a critical value \t\ < t c ~ 1, the dispersion-minimum of the lowest band lies 
at the edges of the Brillouin zone, i.e. at q = ±Qi = (0, ijr^) fc> r our gauge choice with 
a hexagonal 12-site unit cell. In this case eight modes become critical at the confinement 
transition and the resulting unit cell has 36 sites. The corresponding eigenvalue of the 
interaction matrix J±q ± is VQ and the eight eigenvectors occur in complex conjugate pairs 
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Qi 
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Note that in this case the magnetization at lattice site R = 2nu + 2mv and sublattice site 
j is given by 

0,(R) = e^ R ^l + c.c. (3.9) 



n=l..A 
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The PSG transformations of the four complex mode amplitudes ip n can be determined 
similarly to the ferromagnetic case. Quite generally they are defined by 



0(/> j (R) = Re\e i <* 1 < 6K ) ^6v%l] = Re[e lQlR £ (<^) n ug 



71=1. ..4 



71=1. ..4 



(3.10) 



If we define the vector \1/ = (^i, .., tp^, .., tplj , the PSG transformation matrices corre- 
sponding to translations by u, reflections about the x-axis and rotations by n/3 around the 
central site, that act on the vector ^ take the form 
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(3.11) 



(3.12) 



(3.13) 



As in the previous subsection, we used the GAP program to determine that these three 
matrices generate a 288 element subgroup of 0(8) which is isomorphic to GL(2,Z 3 ) x D 3 
The Ginzburg-Landau functional is again given by all homogeneous polynomials that are 
invariant under this group. At fourth order there are five such polynomials, thus the GL- 
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functional depends on the coupling constants r,u, ai, a 4 and is given by 
£4 = H£ + u ^n) + [«i - 2a 4 cos(2(0 2 - 3 )) 

n=l.A 



1 - 20 2 + 4 ) - V3a 2 sinful — 9a) 
VSa 2 sin(0i — 4 )J 



ai 

— a 2 cos(0i — 4 ) — 2a 3 cos 
a 2 cos(6»i - 64) + 2a 3 cos(0i - 20 3 + 4 ) + 
ai + a 4 cos(2(^ - 4 )) - v 7 ^ sin(2(0i - 4 )) 

a 2 cos(0i — 2 ) + 2a 3 cos(0i + 6 2 — 20 4 ) — V3a 2 sin(0! — 6 2 
+ a 4 cos(2(#! - 6 2 )) + VSa 4 sin(2(0i - 2 )) 



ai 



2a 3 cos(20x — 2 — 64) + a 2 cos(0 2 — 64) — v / 3a 2 sin(0 2 — 64) 
ai + a 4 cos(2(0 2 - 4 )) + sin(2(0 2 - 4 )) 

- 20 3 + 6 4 ) + v / 3a 2 sin(0 2 - 4 ) 



+ip2^lip4 - a 2 cos(0 2 - 9 A ) + a 3 cos(0 2 

+V3a 3 sin(0 2 - 20 3 + 64) 
+ipiip 2 ipl - a 2 cos(0i - 2 ) + «3 cos(0i + 2 - 20 3 ) + A/3a 2 sin(0i - 2 ) 

-«s/3a 3 sin(0! + 2 - 20 3 ) 
+M3^l a 2 cos(0i - 3 ) + 2a 3 cos(0i + 3 - 20 4 ) - VSa 2 sin(0i - 
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- a 3 cos (20 2 - 3 - 4 ) + a 2 
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2 



- 2a 3 cos(20i - 3 - 4 ) - a 2 cos(0 3 - 4 ) + V3a 2 sin(0 3 - 4 )] 
ai + a 4 cos(2(0 3 - 4 )) + y/3a 4 sin(2(0 3 - 4 )) 

- a 2 cos(0i - 3 ) + a 3 cos(0i - 20 2 + 3 ) + VSa 2 sin(0i - 3 ) 
-VSa 3 sin(0! - 20 2 + 3 ) 

- 2a 2 cos(0 2 - 3 ) - a 3 cos(0 2 + 3 ) - V3a 3 sin(0 2 + 3 ) 



+^1^2^3 



+^2^3^. 



(3.14) 



Here we have expressed the complex mode amplitudes in terms of their absolute value and 
phase ip n e ten . Note that this fourth order Ginzburg-Landau functional has a remaining 
continuous U(l) symmetry, since it is invariant under a change of all phases n — > 9 n + x- 
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FIG. 11: Phase diagram obtained from the GL-functional (3.15) as a function of the couplings ci2 
and 04. The other parameters are fixed at u = 1, a\ = 1/2, 03 = 02, «5 = 1/20 and ciq = —1/25. 
The different phases are described in the text. 

The "magnetization" (3.9) is not invariant under this U(l) transformation, however. In 
order to break this continuous degeneracy we need to include higher order terms in the 
Ginzburg-Landau functional. Among the invariant sixth order polynomials there are five 
which break the U(l) symmetry. A full analysis of the invariant GL-functional at sixth order 
is beyond the scope of this paper, and we restrict ourselves to the simplest [/(l)-breaking 
sixth order term. The [/(l)-invariant sixth order terms are not expected to qualitatively 
change our results at the confinement transition. 

In the remainder of this section we are going to consider the following fourth order GL- 
functional, including the simplest invariant sixth order [/(l)-breaking polynomial, which 
takes the form 



In total our simplified GL-functional thus has seven coupling constants. Again, a complete 
analysis of the phase diagram as a function of these seven couplings is hardly feasible. A 
representative slice of the phase diagram is shown in Fig. 11, where we have fixed the 
values 43 u = 1, a x = 1/2, = 0-2, 0.5 = 1/20 and = —1/25 and show the different 
phases as function of the two remaining parameters 02 and 04. All phases have a 36-site 
unit cell and are invariant with respect to translations by 4u — 2v and 4v — 2u. In the 
above mentioned parameter regime there are four different ordered phases. VBS 1a has a 
7r/3 rotational symmetry, but is not reflection symmetric. A dimer representation of this 
state is shown in Fig. 1, which was obtained by putting a dimer on every kagome bond 
that intersects a frustrated bond on the dice lattice. This dimer covering suggests that 
our VBS 1a state is identical to previously found valence bond solid states on the kagome 
lattice which maximize the number of perfectly flippable hexagons 4,27 . The VBS 2a phase 




(3.15) 



n=1...4 
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FIG. 12: Bond patterns in the VBS 1a, 2a, 3a and 4a phase (clockwise, starting from the upper 
left). Again we plot the gauge invariant bond-strength Jij(f)i<j)j for nearest neighbor bonds on the 
dice lattice, mapped to the corresponding kagome bonds. Black dashed lines represent satisfied 
bonds {—Jij(j>i(j>j < 0), red lines are frustrated bonds. The thickness of the lines is proportional to 
the bond strength. 

is symmetric under 2ir/3 rotations and has a reflection symmetry. No rotational symmetry 
is present in the VBS 3a and 4a phases. The 3a phase has a reflection symmetry, however, 
which is not present in the 4a phase. Bond patterns of all four phases are shown in Fig. 12. 

We performed a one- loop RG calculation for (3.14) and found six additional fixed points 
besides the Gaussian one, but all of them turned out to be unstable. Again it would be 
useful to revisit this issue using higher loop methods. 37 

IV. CONCLUSIONS 

We have studied confinement transitions of Z2 spin liquids of Heisenberg antiferromagnets 
on the kagome lattice by constructing field theories that are consistent with the projective 
symmetry group of the vison excitations. Depending on the sign of the next-nearest neigh- 
bor interaction between the visons, we found that the visons transformed under the group 
GL(2,Z 3 ) for the simpler case, and under GL(2,Z 3 ) x D 3 for the other case. Our analysis 
shows that possible VBS phases close to the confinement transition are strongly constrained 
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by the vison PSG. We found VBS states that break the translational symmetry of the 
kagome lattice with either 12- or 36-site unit cells, for the two vison PSGs respectively 
The two possible VBS states with 12-site unit cells do not break the rotation symmetry 
of the kagome lattice but one of them breaks the reflection symmetry; this state is closely 
connected to the "diamond pattern" enhancement observed in the recent numerical study of 
Yan et al. 8 As far as possible VBS states with 36-site unit cells are concerned, our analysis is 
not exhaustive. Nevertheless, we found different VBS states with either full, reduced or no 
rotation symmetry, as well as states that do or do not break the reflection symmetry of the 
kagome lattice. Our results should be useful in more completely characterizing spin liquids 
in numerical or experimental studies of the kagome antiferromagnet. 

Analogous analyses for Z 2 spin liquids on other lattices have been carried out in other 
cases (see Appendix B). In all other cases, the effective theory for confining transition has 
an emergent continuous symmetry, and the criticality can be computed using properties of 
the Wilson-Fisher fixed point; reduction to the discrete lattice symmetry appears only upon 
including higher-order couplings which are formally "irrelevant" at the critical fixed point. 
The kagome lattice is therefore the unique case (so far) in which the reduction to discrete 
lattice symmetry appears already in the critical theory: these is the theory in Eqs. (3.8), and 
its relevant quartic couplings are invariant only under discrete symmetries. This suggests 
that the numerical studies of confinement transitions may be easier on the kagome lattice. 
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Appendix A: Berry phase of a vison 

This appendix will compute the Berry phase of a vison moving around a S — 1/2 spin of 
the antiferromagnet, as illustrated in Fig. 4. We will follow the method of Section III. A of 
Ref. 31, generalized to a Z 2 spin liquid as in Ref. 11. 

We consider the time-dependent Schwinger boson Hamiltonian 

ni(r) = -J2 QW e °f> b t* b if> + H - c - + E K(r)b\ a b ia , (Al) 

i<j i 

where the r dependence of and is chosen so that the vison executes the motion shown 
in Fig. 4, while always maintaining the constraint in Eq. (2.2). 
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We compute the Berry phase by working with the instantaneous ground state of %&(t). 
This is facilitated by a diagonalization of the Hamiltonian by performing a Bogoliubov 
transformation to a set of canonical Bose operators, 7^ a , where the index [i = 1 . . . N s , 
where N s is the number of lattice sites. These are related to the b ia by 

h ia = {U^(r) lm - V^r)e a ^) . (A2) 

The N s x N s matrices £/^(r), Uj M (r) perform the Bogoliubov transformation, and obey the 
following identities: 11 





= 1 
= 1 
= 

UV* + V*U T = 0, (A3) 

where a) is a diagonal matrix containing the excitation energies of the Bogoliubov quasipar- 
ticles, and all quantities in Eq. (A3) have an implicit r dependence. 

We can use the above transformations to write down the instantaneous (unnormalized) 
wavefunction of the vison as the unique state which obeys \ ^/ v ) = for all /i, a: 

m = zw{^n 3 e aP b\ a b))j |0>, (A4) 

where 

/5 = E ( u ~ l X ^)„, ■ < A5 > 

Then the Berry phase accumulated during the r variation of Hb( T ) is 



Im(^| |^) = ilmTr 



x 1 dr 



\ dr dr 



(A6) 



We now assume that the Hamiltonian in Eq. (Al) preserves time-reversal symmetry. 
Then, we can always choose a gauge in which the Qij, U ipi and are all real. Under these 
conditions, the expression in Eq. (A6) vanishes identically. It is clear that this argument 
generalizes to the case where we project the wavefunction to boson states which obey the 
constraint in Eq. (2.2). 

We have now shown that no instantaneous Berry phase is accumulated during the vison 
motion of Fig. 4. Under these conditions, the total gauge-invariant Berry phase is simply 
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FIG. 13: Confining phase on the triangular lattice. Plotted is the gauge invariant bond-strength 
Jij4>i4>j f° r nearest neighbor bonds on the frustrated honeycomb lattice, shown on the corresponding 
triangular lattice bonds. Black lines represent satisfied bonds (—Jij4>i(f)j < 0). There are no 
frustrated bonds in the confining phase. The thickness of the bonds is proportional to the bond- 
strength. The two different patterns arise due to a crossover when the sign of the 0(4)-breaking 
term in the GL-functional changes. 

equal to the phase difference between the wavefunctions in the initial and final states. 9 As 
shown in Fig. 4, this phase difference is it. 

Appendix B: Effective Ising models for visons on various other lattice geometries 

In this Appendix we summarize the results of a Ginzburg-Landau analysis of frustrated 
transverse field Ising models (TIMs) on various lattice geometries. These models describe the 
low energy properties of different frustrated Heisenberg antiferromagnets in terms of their 
vison exciations. Some of these results have been discussed previously in the literature. 

1. Triangular lattice 

The vison excitations of a Heisenberg antiferromagnet on the triangular lattice are de- 
scribed by a frustrated TIM on the dual honeycomb lattice, which has been studied previ- 
ously by Moessner and Sondhi 25 . A Ginzburg-Landau analysis reveals four critical modes 
and the corresponding PSG transformation matrices generate a 288 element subgroup of 
0(4) which is isomorphic to (C3 x GL(2, Z3)) x C2, 33 where C n denotes the cyclic group of 
order n. PSG matrices for a specific gauge choice can be found in Ref. 25. An 0(4)-breaking 
term appears at sixth order in the GL-functional, the minimization of which gives rise to 
a single confined phase with a 24-site unit cell (i.e. a 12-site unit cell on the triangular 
lattice) that is symmetric under 27r/3-rotations and reflections. Bond patterns of this phase 
are shown in Fig. 13. Note that there is a transition when the sign of the 0(4)-breaking 
term changes. 
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2. Honeycomb lattice 

For the frustrated honeycomb lattice antiferromagnet the visons are described by an 
antiferromagnetic TIM on the dual triangular lattice. There are two critical modes at the 
Brillouin zone edges Q = ±(4-7r/3, 0) and the PSG matrices corresponding to translations 
T by any basis vector, rotations Rq and reflections I y about the y-axis of the two mode 
amplitudes are given by 

T = -(l 2 + iV3a z )/2 , R 6 = I y = C r x (Bl) 

where Oi denote the Pauli matrices. These PSG matrices generate the 6 element dihedral 
group D%, i.e. the symmetry group of the equilateral triangle. The invariant Ginzburg- 
Landau functional has been discussed previously by Blankenschtein et al. M , who showed that 
an 0(2) symmetry breaking term appears at sixth order. 23 Minimizing the GL functional 
gives rise to only one possible confining phase which breaks the translational symmetry. For 
a particular sign of the sixth order term, 23 the confining phase has a three-site unit cell 
(i.e. six sites per unit cell on the honeycomb lattice) and is symmetric both with respect to 
rotations and reflections; the corresponding dimer pattern on the honeycomb lattice has the 
maximal number of one perfectly flippable hexagon per six-site unit cell 25 and is identical 
to the VBS state found in Ref. 31. A plaquette-like phase is obtained for the other sign of 
the sixth-order term. 23 More complex minima structure for the vison dispersion have also 
been considered in Ref. 23. 

3. Square lattice 

The effective vison model for the frustrated square lattice Heisenberg antiferromagnet is 
a frustrated TIM on the dual square lattice. 21,45 In an appropriate gauge the two critical 
modes appear at zero momentum 46 and the corresponding (gauge-dependent) PSG-matrices 
for translations T x , rotations i?4 and reflections I x take the form 

T x =(a x + a z )/V2, Rt = a z , h = 1 2 ■ (B2) 

These matrices generate the 16 element dihedral group D 8 . An invariant GL-polynomial 
that breaks the 0(2) symmetry appears only at eight order, as discussed by Blankenschtein 
et al. A1 . Depending on the sign of this eight order term, two different confining phases are 
possible. Both phases are reflection symmetric, break the translational symmetries and have 
a four-site unit cell. One of the two phases is invariant under 7r/2-rotations, whereas the other 
one has a reduced rotational symmetry and is only invariant with respect to 7r-rotations. 
These are the familiar 'plaquette' and 'columnar' VBS states. 48 

More complex vison dispersion structures, with further-neighbor couplings, have been 
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described recently in Ref. 23. 
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